Combined analysis of rj meson hadro- and photo-production off nucleons 



00 

o 
o 

(N 

Or 
< 
(N 



43 

Or 

43 



> 

On 

rn 
cn 

o 

oo 
O 



X 



K. Nakayama, 1,2 '* Yongseok Oh, 3, t and H. Haberzettl 4, * 

1 Department of Physics and Astronomy, University of Georgia, Athens, GA 30602, U.S.A. 
2 Institut fur Kernphysik (Theorie), Forschungszentrum Julich, 52425 Jiilich, Germany 
''Cyclotron Institute, Texas A&M University, College Station, TX 77843, U.S.A. 
* Center for Nuclear Studies, Department of Physics, The George Washington University, Washington, DC 20052, 

(Dated: 30 March 2008) 

The r/-meson production in photon- and hadron-induced reactions, namely, 7p — » pr), ir~p —* nr\, 
pp — ► ppn, and pn — > pun, are investigated in a combined analysis in order to learn about the relevant 
production mechanisms and the possible role of nucleon resonances in these reactions. We consider 
the nucleonic, mesonic, and nucleon resonance currents constructed within an effective Lagrangian 
approach and compare the results with the available data for cross sections and spin asymmetries 
for these reactions. We found that the reaction 7p — * pr] can be described well with the inclusion 
of the well-established Su(1535), S , n(1650), Di3(1520), and Di3(1700) resonances, in addition to 
the mesonic current. Consideration of other well-established resonances in the same mass region, 
including the spin-5/2 resonances, Dis(1675) and Fis(1680), does not further improve the results 
qualitatively. For the reaction ir~p — > nrj, the Pi3(1720) resonance is found to be important for 
reproducing the structure observed in the differential cross section data. Our model also improves 
the description of the reaction NN — > NNrj to a large extent compared to the earlier results by 
Nakayama et al. [Phys. Rev. C 68, 045201 (2003)]. Further improvement in the description of these 
reactions and the difficulty to uniquely determine the nucleon resonance parameters in the present 
type of analysis are discussed. 

PACS numbers: 25.20.Lj, 13.60.Le, 13.75.Gx, 14.20.Gk 
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I. INTRODUCTION 

One of the primary motivations for studying the pro- 
duction of mesons off nucleons is to investigate the struc- 
ture and properties of nucleon resonances and, in the case 
of heavy-meson productions, to learn about hadron dy- 
namics at short range. In particular, a clear understand- 
ing of the production mechanisms of mesons heavier than 
the pion still requires further theoretical and experimen- 
tal investigation. Apart from pion production, the ma- 
jority of theoretical investigations of meson-production 
processes are performed within phcnomcnological meson- 
exchange approaches. Such an approach forces us to cor- 
relate as many independent processes as possible within 
a single model if one wishes to extract meaningful physics 
information. Indeed, this is the basic motivation behind 
the coupled-channels approaches. 

In this paper, we present the result of our investigation 
of ry-meson production in both the photon- and hadron- 
induced reactions, 
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More specifically, we perform a combined analysis of the 
reactions, jp — > pr/, rr~p — > nrj, pp — + pprj, and pn — ► pnn 
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based on an effective Lagrangian approach. 

The amount of data available for ?7-meson production 
is now considerable. In particular, in photoproduction 
processes off protons, a new generation of high-precision 
data are now available, not only for total and differen- 
tial cross sections in a wide range of energies starting 
from threshold [1] but also for beam and target asym- 
metries [2], Taken in combination, these data sets offer 
better opportunities for investigating the properties of 
nucleon resonances. In particular, much more detailed 
studies than in the past are possible now for resonances 
that may perhaps couple strongly to Nrj but only weakly 
to Ntt. In this respect, the recent data on the quasi- 
free jn — > nrj process [3] have attracted much interest in 
77-production processes in connection to the possible exis- 
tence of a narrow (crypto-exotic) baryon resonance with 
a mass near 1.68 GeV which is still under debate [4]. To 
date, the narrow bump structure in the jn — > nq reaction 
cross section was reported by several groups [5, 6]. The 
existence of a narrow baryon resonance near 1.68 GeV 
has been also suggested earlier by several authors [7]. 

There are a large number of theoretical investigations 
of r) photoproduction, mostly from the early 1990s to the 
present [8] and, especially, off protons. Most of them fo- 
cus on the role of nucleon resonances and the extraction 
of the corresponding resonance parameters, but a vari- 
ety of issues have also been addressed, such as the NNrj 
coupling constant [9], the UaO-) anomaly [10], and the 
extended chiral symmetry [11]. Among the more recent 
calculations which analyze the recent high-precision data 
arc those of Refs. [12-14]. The 77-MAID approach [12] 
is an isobar model which includes, in addition to the t- 
channel vector-meson exchange or its Reggeized version, 
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a set of well-established spin- 1/2, -3/2, and -5/2 reso- 
nances. This model has been applied to the analyses of 
data for photoproduction as well as for electroproduction 
of the r\ meson. The Bonn-Gatchina-Giessen group [13] 
has developed a model which has been employed in a 
combined partial wave analysis of the photoproduction 
data with Nit, Nr], AK, and T,K final states. This model 
considers fourteen nuclcon resonances and seven Delta 
resonances for achieving a reasonable overall agreement 
with the whole database considered in their analysis. 1 
He, Saghai, and Li have analyzed the 77-photoproduction 
reaction in a chiral constituent quark model [14] by con- 
sidering all of the one- to four-star rated resonances listed 
in the Review of Particle Data Group (PDG) [15]. 

In hadronic reactions, a noticeable amount of data for 
NN — ► NNi] has been accumulated [16-19]. Here we 
have pp and pi] invariant mass distributions and the an- 
alyzing power near threshold in addition to differential 
and total cross sections. The NN — > NNr] process is 
particularly relevant for studying the role of the Nr] final- 
state interaction (FSI). Most of the existing calculations 
(see, e.g., Refs. [20, 21] and references therein; see also 
Ref. [22]) take into account the effects of the NN FSI 
in one way or another which is well-known to influence 
the energy dependence of the cross section near thresh- 
old. Calculations which include the Nrj FSI to lowest 
order [23] reproduce the bulk of the energy dependence 
exhibited by the data. However, they are not sufficient to 
reproduce the pp invariant mass distribution measured by 
the COSY-TOF [17] and COSY-11 Collaborations [18]. 
Thus, in order to explain the observed pp invariant mass 
distribution, the importance of the three-body nature of 
the final state (in the S'-wave) has been emphasized [24] 
or an extra energy dependence in the basic production 
amplitude has been suggested [25]. Yet, another possi- 
bility has been offered, which is based on the higher par- 
tial wave (P-wave) contribution [26]. We observe that 
what is actually required to reproduce the measured pp 
invariant mass distribution is an extra p' 2 dependence, 
where p' denotes the relative momentum of the final pp 
subsystem. Obviously, this can be achieved either by an 
S-wave or by a P-wave contribution. Note that the NN 
P-wave ( 3 Po) can also yield a flat proton angular distri- 
bution as observed in the corresponding data. The model 
of Ref. [26], however, underpredicts the measured total 
cross near threshold to a large extent. One of the objec- 
tives of the present study is to resolve this discrepancy. 
In any case, as pointed out in Ref. [26] , the measurement 
of the spin-correlation functions should help settle the 
question of the S-wave versus P-wave contributions in a 
model-independent way. 

Most data for more basic (two-body) irN — * Nr/ re- 
action have been obtained in the 1960s through 1970s; 



The A resonances do not coupled to the Nr) and AK channels 
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they are rather scarce and less accurate [27] than for the 
other two reactions mentioned above. Recently, the Crys- 
tal Ball Collaboration has measured the differential and 
total cross sections of this reaction near threshold [28]. 
Theoretically, this reaction has been studied mostly in 
conjunction with other reactions in a combined analy- 
sis [29] or in coupled-channels approaches [30-34] in or- 
der to constrain some of the model parameters. Recently, 
Zhong et al. [35] have extended the chiral constituent 
quark model for meson photoproduction [14] to this re- 
action. Also, irN — > Nrj has been studied within a heavy- 
baryon chiral perturbation theory [36]. Arndt et al. [37] 
have investigated the role of irN — > Nr] on the Nr] scat- 
tering length within a coupled-channels analysis of this 
reaction and the elastic ttN scattering. 

In the present work, we consider the three reactions 
mentioned above in the following manner: The photo- 
production reaction is calculated by considering the s-, 
u- and t-channel Feynman diagrams plus a generalized 
contact term [38] , which ensures gauge invariance of the 
total amplitude, in addition to accounting for the final- 
state interaction effects. (See Ref. [39] for details.) The 
ttN — > Nr] reaction is calculated in the tree-level ap- 
proximation including the s-, u-, and i-channels. To the 
extent that this reaction is dominated by the excitation 
of the Sii(1535) resonance at least for energies close to 
the threshold, this should be a reasonable approxima- 
tion if we confine ourselves to energies not too far from 
the threshold. The NN — > NNr] process is calculated 
in the Distorted- Wave Born Approximation (DWBA), 
where both the NN FSI and the initial-state interaction 
(ISI) are taken into account explicitly [20]. The NN FSI 
is known to be responsible for the dominant energy de- 
pendence observed in the total cross section (apart from 
the dependence due to the phase space) arising from the 
very strong interaction in the S'-wave states at very low 
energies [43]. As for the basic meson-production am- 
plitude, our model includes the nucleonic, mesonic, and 
nucleon resonance currents which are derived from the 
relevant effective Lagrangians. 

Ultimately, our goal is to perform a more complete 
model calculation where the relevant FSIs are taken into 
account explicitly. However, before being able to un- 
dertake a complex calculation that couples many chan- 
nels, we need to learn some of the basic features of me- 
son production (in particular, those of ry-meson produc- 
tion) within a simplified model where these basic features 
may be revealed and analyzed in a much easier man- 
ner. In this regard, one of the major purposes of the 
present investigation is to show that consideration of the 
hadronic reactions NN — > NNM , in conjunction with 
more basic two-body reactions, would greatly help in the 
study of nucleon resonances, especially, in imposing much 
stricter constraints on the extracted resonance-nucleon- 
meson (RNM) coupling strength involving a meson M 
other than the pion. In fact, currently, our knowledge of 
the branching ratios of the majority of the known reso- 
nances is very limited [15]. 
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This paper is organized as follows. In the next section, 
we briefly describe our model for the reactions listed in 
Eq. (1). The results of the corresponding model calcula- 
tions are presented and discussed in Sec. III. Section IV 
contains a summary and conclusion. Some details of the 
present model are given in the Appendix. 



II. REACTION MECHANISMS 

In the present work, the ?7-meson production pro- 
cesses are treated within a relativistic meson-exchange 
approach, whose dynamical contents are summarized by 
the Feynman diagrams displayed in Figs. 1,2, and 3, for 
the reactions jN -> Nrj, ttN -> Nr], and NN -> NNrj, 
respectively. We employ phenomenological form factors 
at hadronic vertices to account for the structure of the 
corresponding hadrons. 

For the photoproduction process, the total amplitude 
in the present work is given by the Feynman diagrams 
displayed in Fig. 1. Gauge invariance of the total am- 
plitude is ensured by a generalized contact current given 
in Refs. [38, 39] that follows the general formalism of 
Refs. [40-42]. This contact term provides a rough phe- 
nomenological description of the FSI which is not treated 
explicitly here [39]. The details of the present approach 
are fully described in Ref. [38], where r/-meson produc- 
tion in photon- and hadron-induced reactions was inves- 
tigated, and they will not be repeated here. One new 
feature in the present work, however, is the inclusion of 
spin-5/2 resonance contributions. 



As mentioned in Sec. I, the reaction irN — > Nr) is calcu- 
lated within the tree-level approximation. To the extent 
that it is dominated by the Sii(1535) resonance contri- 
bution, this is a reasonable approximation at least near 
threshold. The total amplitude for this reaction is, there- 
fore, given by the Feynman diagrams displayed in Fig. 2. 

As for the NN — > NNr) reaction, the total amplitude 
is calculated within the DWBA, 

M =(T f G f + l)J(l + G j T i ), (2) 

where T i(/) denote the NN ISI (FSI); G l{f) stands for the 
corresponding JVJV propagator; J denotes the basic pro- 
duction amplitude displayed in Fig. 3 and is constructed 
from the interaction Lagrangians given in the Appendix. 




FIG. 2: Feynman diagrams contributing to 7r~p — * nr). The 
notation is the same as in Fig. 1. The nucleonic current (nuc) 
referred to in the text corresponds to the first two diagrams 
on the top line; the meson-exchange current (mec) and res- 
onance current contributions correspond, respectively, to the 
rightmost diagram on the top line and the two diagrams on 
the bottom line. 




FIG. 1: Feynman diagrams contributing to jp — ► prj. Time 
proceeds from right to left. The wavy, solid, and dashed lines 
represent the photon, the nucleon, and the r) meson, respec- 
tively. The intermediate baryon states are denoted by N and 
R for the nucleon and the nucleon resonances, respectively. 
The intermediate mesons in the f-channel include the p and 
u>. The external legs are labeled by the four-momenta of the 
respective particles and the labels s, u, and t of the hadronic 
vertices correspond to the off-shell Mandelstam variables of 
the respective intermediate particles. The three diagrams in 
the lower part are transverse individually; the three diagrams 
in the upper part are made gauge-invariant by the general- 
ized contact current depicted in the top-right diagram. The 
nucleonic current (nuc) referred to in the text corresponds to 
the top line of diagrams; the meson-exchange current (mec) 
and resonance current contributions correspond, respectively, 
to the leftmost diagram and the two diagrams on the right of 
the bottom line of diagrams. 




+ R |;M + Mil R 

pi 6 — Pi p'i — 6 Pi 




FIG. 3: Basic production amplitude for NN -> NNrj. Time 
proceeds from right to left. As in Fig. 1, N and R denote 
the intermediate nucleon and resonances, respectively, and 
M incorporates all exchanges of mesons n, r), p, u>, a, and 
do (former <5) for the nucleon graphs and tv, n, p, and w for 
the resonance graphs. External legs are labeled by the four- 
momenta of the respective particles; the hadronic vertices s, 
u, and t correspond to the same kinematic situations, respec- 
tively, as those identified similarly in Fig. 1. The nucleonic, 
resonance, and meson-exchange contributions referred to in 
the text correspond, respectively, to the first, second, and 
third lines of the Feynman diagrams on the right-hand side. 
Diagrams with p 1 «-» p 2 are not displayed but are included 
in the calculation. The full amplitude, including the NN ISI 
and FSI contributions, is given by Eq. (2). 
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Further details of the present approach to this reaction 
can be found in Ref. [20]. Also, we use the NN interac- 
tion based on the Paris potential [44] , which includes the 
Coulomb interaction as well [45]. 

In the Appendix, we present all the hadronic and 
electromagnetic interaction Lagrangians and propaga- 
tors necessary for computing the diagrams displayed in 
Figs. 1,2, and 3 within the present approach. The phe- 
nomenological form factors used in this model are also 
given in the Appendix. 

The free parameters of our model — the resonance 
parameters, the NNr] coupling constant, and the cut- 
off parameter A* at the electromagnetic vector-meson 
exchange vertex — are fixed such as to reproduce the 
available data in a global fitting procedure of the three 
reaction processes listed in Eq. (1). 



III. RESULTS AND DISCUSSION 

In this section, we present and discuss the results of our 
model calculation. The basic strategy of our approach is 
in principle the same as that of Ref. [38] , namely, we start 
with the nuclcon plus meson-exchange currents and add 
resonance contributions one by one as needed in the fit- 
ting procedure until achieving a reasonable description of 
the available experimental data for the reactions listed in 
Eq. (1). Apart from the dominant 5n(1535) resonance, 
we allow for other well-established resonances of spin-1/2, 
-3/2, and -5/2 in this model. We find that, in addition 
to spin-1/2 resonances, one needs at least spin-3/2 res- 
onances — in particular, D\% resonances — in order to 
obtain a reasonable description of the data in the present 
approach. 

Following Ref. [38], for each resonance we take into 
account only the branching ratios P Nw , and P N7T7T 

corresponding to the respective hadronic decay channels 
R — ► Ntt, R — > Nt], and R — > Nnn. The latter accounts 
effectively for all the other open decay channels. Note 
that the branching ratios /3 Nn and /3 N are related to 
the corresponding RNn and RNr] coupling constants in 
the interaction Lagrangians and, as such, they are not 
free parameters of the model. The same is true for the 
branching ratio "f Nl , associated with the radiative decay 
channel, which is related to the corresponding RNj cou- 
pling constants. Then, in view of the constraint given 
by Eq. (A. 31), the branching ratio /3 N7T7T is not a free 
parameter, either. Here, we emphasize that, in contrast 
to Ref. [38] — where some assumptions were made con- 
cerning the values of the branching ratios (3 Nn — no such 
assumptions are enforced in the present work since the si- 
multaneous consideration of the reaction processes listed 
in Eq. (1) allows us, in principle, to extract the i?A^7r and 
RNr/ coupling constants separately. 

The coupling constants of the RNV (V = p,ui) in- 
teraction Lagrangians are required in the calculation of 
the NN — > NNr/ reaction. Therefore, in principle, one 
should also take into account the hadronic branching ra- 



tios involving vector mesons (such as the p and u>). How- 
ever, since we have restricted ourselves to nucleon reso- 
nances for which the decay channels R ~ * NV are either 
closed or nearly closed, we have set the associated branch- 
ing ratios to zero (/3 NV — 0). 2 In a test calculation, we 
found that this does not affect our conclusions. Obvi- 
ously, in a more refined calculation, this condition should 
be relaxed. In a more complete (coupled-channels) dy- 
namical model approach, the branching ratios and total 
widths will be generated by the model via the dressing 
mechanism of the corresponding vertices and resonance 
masses. 

The results shown here are not necessarily of the best 
fit quality one can achieve within the present approach. 
Rather, they are sample fit results that illustrate different 
dynamical features one may obtain in this type of anal- 
ysis. The resonance parameters are obtained by global 
fitting to the available data for the reactions mentioned 
above. In the following, we present and discuss the re- 
sults for each reaction. 



A. 7p — ► pr] 

We first consider the model in which only the 
5n(1535), 511(1650), £>i 3 (1520), and L>i 3 (1700) reso- 
nance currents are considered in addition to the nuclc- 
onic and mesonic currents (cf. Fig. 1). We find that this 
comprises the minimal set of resonances that are required 
to achieve a reasonable description of the reaction pro- 
cesses listed in Eq. (1). The resulting fitted parameters 
are given in Table I, which are obtained by x 2 fitting 
to the experimental data shown in Fig. 4 (however, not 
taking into account the total cross sections for jp — > pn) . 
Note that for photoproduction, in addition to the electro- 
magnetic couplings, only the hadronic vertices involving 
the r] meson are required for the present calculation. In 
the Tr~p — ► ni] reaction, only the hadronic vertices involv- 
ing i] and 7r are needed, while all the hadronic vertices 
given in the Appendix are required for the calculation of 
the NN — > NNn reaction. 

In Table I, the parameter values in boldface are fixed 
and are not allowed to vary during the fitting procedure. 
The pure pscudovector coupling choice (A = 0) at the 
NNij vertex was motivated by the massless chiral limit 
of the r] meson. (See the Appendix for the definition 
of A.) However, we will also consider the pseudoscalar 
coupling choice (A = 1) in subsection IIIB. Also, the 
most general form of the RNV vertex for spin-3/2 and 
-5/2 resonances involves three independent coupling con- 
stants as exhibited in Eqs. (A. 8) and (A. 10). At present, 
however, information on the corresponding coupling con- 
stants are extremely scarce, especially, on g^ NV defined 



2 Strictly speaking, even these resonances can have non-vanishing 
branching ratios to NV channels due to their large widths. 



■5 



TABLE I: Model parameters fitted to the reactions listed in Eq. (1). (See the Appendix for explanation of the parameters.) 
Values in boldface are not fitted. The branching ratios and /3 N j (j — ir,r),iT-K) are not free parameters but are extracted 
from the corresponding coupling constants, except for f3 N7T7I which is obtained from Eq. (A. 31). The values in square brackets 
are the range estimates quoted in PDG [15]. The data set for jp — ► pr\, ir~p — > nr\, and NN — > NNn are used in the fit. 



Nucleonic current: 










(SiVJVr,! A) 


(0.007, 0.0) 








Mesonic current: 










A; (MeV) 


1168 








An current: 


Sn(1535) 


PDG 


Sn (1650) 


PDG 


M R (MeV) 


1540 


[1525-1545] 


1615 


[1645-1670] 




0.81 




0.52 




(QrNtvi A) 


(-0.67, 0.04) 




(-0.50, 0.77) 




{9RNrii A) 


(-2.61, 0.46) 




(0.98, 0.91) 




(9RNuJ1 /flJVo)) 


(38.52, -1.39) 




(86.78, 6.19) 






(-15.07, 2.67) 




(-39.01, 10.84) 




F R (MeV) 


200 


[125-175] 


144 


[145-185] 


7iv 7 (%) 


0.31 


[0.15-0.35] 


0.24 


[0.04-0.18] 


Pn* (%) 


32 


[35-55] 


27 


[60-95] 


(%) 


65 


[45-60] 


19 


[3-10] 


Avttjt (%) 


2 


[1-10] 


54 


[10-20] 


A13 current: 


Di 3 (1520) 


PDG 


£>i 3 (1700) 


PDG 


M R (MeV) 


1520 


[1515-1525] 


1700 


[1650-1750] 


(9rn~{ ' 9rn~/ ) 


(1.22, -0.18) 




(0.21, -0.22) 






-1.96 




0.58 






-3.32 




-2.40 




fo (1) o (2) o (3) ) 


(-289.94, 242.59, 0.0) 




(341.79, -362.81, 0.0) 




RN pi 9 RN pi 9 RN p) 


(54.58, 154.53, 0.0) 




(-1.61, 90.18, 0.0) 




F R (MeV) 


108 


[100-125] 


94 


[50-150] 


7iv 7 (%) 


0.41 


[0.46-0.56] 


0.07 


[0.01-0.05] 


Av. (%) 


95 


[55-65] 


28 


[5-15] 


/J*„ (%) 


0.03 


[0.23 ± 0.04] 


2 


[0.0 ± 1.0] 


Avtttt (%) 


5 


[40-50] 


70 


[85-95] 



in Eqs. (A. 8) and (A. 10). In the present work, therefore, 
we allow only two structures at those vertices by setting 
the coupling constant g^ NV to zero. 

In this work, we use the resonance masses from the cen- 
troid values quoted in PDG [15]. The exception to this 
is the masses of the 5n(1535) and 511(1650) resonances. 
They were allowed to vary in the fitting process in order 
to reproduce the large cross-section peak exhibited by the 
photoproduction data near the threshold. The quantities 
in the square brackets in Table I are the range estimates 
quoted in PDG [15] and are given here for an easy com- 
parison with the fitted values extracted in the present 
model. As one can see, some parameter values are con- 
siderably outside the range quoted in PDG. However, we 
note that, although various parameters are highly corre- 
lated to each other, the currently available data used in 
the present study are not sufficient to uniquely constrain 
the model parameters. As a result, different parameter 
sets may provide fits of comparable quality; some of them 
are discussed in this paper. 



Figure 4 displays the results for photoproduction ob- 
servablcs corresponding to the parameter set of Table I. 
The total cross section shown in Fig. 4(a) provides the 
line styles used in all four parts. We start with the dis- 
cussion of the rj angular distribution shown in Fig. 4(b). 
As one can see, the flat angular distribution near the 
threshold is dominated by the Su resonance. When 
the photon incident energy (in the laboratory frame) T 7 
is larger than about 1.0 GeV, both the D i3 resonance 
and the mesonic currents become relevant for reproduc- 
ing the shape of the measured differential cross sections. 
As energy increases, the shape becomes more and more 
forward-peaked, which is a well-known feature of the t- 
channcl mesonic current contribution. In the energy re- 
gion of T 7 w 1.1-2.0 GeV, some details of the measured 
angular distribution are still not well explained, which 
indicates that our model should be improved by a more 
refined and quantitative calculation. We leave such an 
investigation to a future work. The nucleonic current 
contribution is negligible because of very small coupling 
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cos (%) cos(e v ) 

(b) (d) 

FIG. 4: (Color online) Results for the reaction 7p — -> pr/ employing the parameters of Table I. (a) Total cross section as a 
function of the total energy of the system W. The line styles identified here apply to all four parts of this figure, (b) r\ angular 
distribution in the center-of-mass frame, (c) Beam asymmetry E_b and (d) target asymmetry Et in the center-of-mass frame. 
The numbers in (b,c,d) are the incident photon laboratory energies T 7 in MeV. In (c,d) only the total results are shown. The 
data are from Refs. [1, 2]. 



constant g^^ n resulting from the fit, which is preferred angles and higher energies. The nucleonic current con- 
by the small angular distribution measured at backward tributes mostly at backward angles and at high energies 
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TABLE II: Another set of fitted parameters. See the caption of Table I for details. 



Nucleonic current: 












(U.Uoo, U.UJ 








Mesonic current: 










A v (MeV) 


1202 








Nn current: 


Sn(1535) 


PDG 


Sn (1650) 


PDG 


Tl /T f TV T T T\ 

M R (MeV) 


1527 


[1525-1545] 


1625 


[1645-1670] 


(i) 
9rn~/ 


—0.55 




0.38 






(0.43, 0.005) 




(-0.37, 0.84) 






(1.92, 0.44) 




(0.68, 0.86) 




{9rn^ i fnNu ) 


(-3.86, 0.79) 




(8.81, -29.85) 




(Si?jvpi /fljvp) 


(29.18, 8.40) 




(4.09, 3.63) 




r fl (MeV) 


99 


[125-175] 


162 


[145-185] 


7iv 7 (%) 


0.27 


[0.15-0.35] 


0.12 


[0.04-0.18] 


1 N ir [JO) 


2b 


[oo-ooj 


13 


[oU-yoj 


PNr/ (/Oj 


63 


[4o-dUJ 


n 
9 


To i nl 

[3-10] 


toy \ 

PNirn 1/0 J 


10 
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through the u-channel diagram (cf. Fig. 1). We will come 
back to this issue later in subsection IIIB. 

Shown in Fig. 4(a) are the total cross section for r\ 
photoproduction obtained with the resonance parameters 
given in Table I. Here we mention that the total cross 
section data for photoproduction were not included in 
the global fitting process. As one can see, the large peak 
that rises sharply from threshold is due to the dominat- 
ing 5ii(1535) and 511(1650) resonances. Although both 
the D13 resonance and mesonic currents are small, their 
contribution beyond W ~ 1.6 GeV cannot be ignored be- 
cause of the interference with the large Sn current con- 
tribution. Around W = 2.0 GeV, all the resonance and 
mesonic currents become comparable to each other, but 
at higher energies, the mesonic current yields the largest 
contribution and the Su resonance contribution becomes 
negligible as expected from their relative low mass. The 
structure shown by the data at W <~ 2.2 GeV seems 
to suggest possible contributions from other higher mass 
resonances. Analyses of this structure and other details 
are left to a future work. 

Figure 4 also shows the results for the beam and tar- 
get asymmetries. We find that, by and large, the beam 



asymmetry (£.b) and the target asymmetry (St) are de- 
scribed reasonably well. 3 The most visible room for im- 
provement exists for the latter, in particular, at lower en- 
ergies, where the data show a sin(20 ?7 ) dependence while 
the model yields nearly flat and vanishing results. This 
is a feature common to all the parameter sets and not 
just to the particular parameter set of Table I. To un- 
derstand the origin of this discrepancy, we consider the 
partial- wave contribution to the spin asymmetries. We 
expect that only a few partial waves contribute mostly 
to this reaction near the threshold. Then, considering 
only the S- and P-waves in the final state and following 
Rcf. [47], the beam and target asymmetries are readily 



3 Quite recently, the GRAAL Collaboration [46] has re-analyzed 
the beam asymmetry in ■yp — > pr) revealing a sharp structure at 
W ~ 1.69 GeV for 8 V ~ 43° — 103°, suggesting the presence of 
a narrow resonance. 
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FIG. 5: (Color online) Same as Fig. 4 corresponding to the fitted parameter set of Table II. 



given by 

^S B = (|a i=1 | 2 -|/? i=1 | 2 )sin 2 (^), 

= Im [6* L=0 (p L=1 ~ a L=1 )} sin(0„) 

+ Im [S* L=1 (p L=1 - a L=1 )] sin(20„), (3) 



which exhibits their respective angular dependence, 
where a L , (3 L , and 5 L stand for the linear combinations 
of the partial-wave matrix elements with the orbital an- 
gular momentum L in the final state. Then, it is easy 
to understand that it is more difficult to describe the 
latter observable because the beam asymmetry does not 
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TABLE III: Same as Table I. The P 13 (1720) resonance is added here to see whether it further improves the fit. No NN -> NNrj 
data were included in this fit. 
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exhibit interference effects while the target asymmetry 
does. Equation (3) also shows that the S'-wave final state 
does not contribute to the beam asymmetry. 

Another parameter set resulting from the global fit, 
employing the same set of nucleon resonances as in Ta- 
ble I, is shown in Table II. This set was obtained by using 
different starting values for the search in parameter space 
during the fit procedure. As a result, many parameter 
values of this set are quite different from those of Ta- 
ble I not only in magnitude but also in relative signs 
for some coupling constants which provides an indica- 
tion of the general reliability of such global fits. The 
corresponding observables are shown in Fig. 5. Here, 
although both the measured total and differential cross 
sections are reproduced with a comparable fit quality to 
the results in Fig. 4, the dynamical content is quite differ- 
ent. In particular, the Z?i3(1700) resonance contribution 
dominates over the other currents in the energy region 
of T 7 w 1.05-1.5 GeV (i.e., W w 1.7-1.9 GeV). It is 
interesting to note that, searching for pole positions of 
the T-matrix in the complex plane as well as performing 
Breit-Wigner parameterizations of resonances, the recent 
partial-wave analysis by Arndt et al. [48] of n ± p elastic 
and charge- exchange processes combined with the reac- 
tion ir~p — > nr\ finds no Z?i3(1700) resonance. Also, some 
of the coupled-channels dynamical models [32, 49] find no 



necessity for this resonance to fit the irN phase shifts. In 
this respect, the present results corresponding to the pa- 
rameter set of Table I — where the contribution of the 
-Di3(1700) resonance is much smaller than that of Ta- 
ble II — are more in line with these findings. It would be 
most interesting to see if the inclusion of the 7^ channel 
in those coupled-channels analyses mentioned above will 
require the L>i3(1700) resonance. The spin observables 
in Fig. 5 exhibit overall the same features as in Fig. 4, 
but for higher energies the target asymmetry shows very 
different angular dependence for the two parameter sets 
which points to the importance of spin asymmetries in 
reducing the ambiguities that otherwise would exist. 

Table III displays a parameter set including the 
Pi3(1720) resonance, in addition to those considered in 
the previous two sets. Here we have considered only the 
IP ~ * PV an d — ► ni] reaction data in the fitting pro- 
cedure. The corresponding results for the observables 
are shown in Fig. 6. As can be seen, the inclusion of 
the Pi3(1720) resonance does not improve significantly 
the description of the data for this photon-induced reac- 
tion. However, as we will see in the following subsection, 
this resonance considerably improves the fit quality of 
the hadronic ir~p — ► nij reaction at higher energies. 

Next, we extend our model by including all of the 
well-established resonances in the mass region of 1500 <~ 
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FIG. 6: (Color online) Same as Fig. 4 corresponding to the parameter set of Table III. 



1700 MeV in order to verify whether these resonances 
can lead to a qualitatively superior description of the data 
compared to the previous cases, where a more limited set 
of resonances has been considered. To this end, we con- 
centrate only on the 72? — > prj reaction. Table IV shows 
the parameter set, where the Z?i5(1675) and Fi5(1680) as 



well as the Pn(1710) resonances are included, in addition 
to those considered in Table III. Here, following Ref. [38], 
we treat the branching ratio f3 NlT as a free parameter to 
be fitted, while the branching ratio f3 N is extracted from 
the product of the coupling constants 9fjNr)9iiN'y m con ~ 
junction with the assumed branching ratio 7^ for the 
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TABLE IV: Same as Table I. More resonances are added here to see whether they would further improve the fit. Here we have 
considered only the 7p — > pr/ reaction in the fitting procedure. 



Nucleonic current: 
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A^is current: 
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radiative decay. The resulting observables with the pa- 
rameter set of Table IV are shown in Fig. 7. We see that 
the overall fit quality does not change significantly from 
the previous sets. Here, the Z?i3(1700) resonance gives 
the largest contribution to the cross section in the energy 
region of W = 1.7-2.0 GeV, a feature similar to what is 
exhibited already in Fig. 5. 



B. NNr/ coupling constant 

In the previous subsection we have shown that the 
present calculations yield very small values of the cou- 
pling constant g NNrj — compatible with zero — due to 
the smallness of the measured cross sections at backward 
angles where T 7 is large. As mentioned before, the 77 
angular distribution becomes very sensitive to g NNr) at 
these kinematics through the u-channel nucleonic current 
contribution. However, one must be cautious in drawing 



conclusions about the extracted value of g NNr) from cal- 
culations based on approaches such as the present one. 
This is due to the fact that we cannot completely dis- 
card the possibility of a relatively large nucleonic current 
contribution interfering destructively with contributions 
from resonances in order to yield 77 angular distributions 
as observed at those kinematics. This point is illustrated 
in Fig. 8, where the results for 77 angular distribution are 
shown together with the data recently reported by Crede 
et al. (CB/ELSA Collaboration) [1] at the two highest 
energies. Figure 8(a) shows the results corresponding 
to the parameter set of Table IV, where the nucleonic 
current contribution is practically zero 4 and cannot be 
seen in the figure. The corresponding nucleon resonance 
contributions are also very small. The mesonic current 



4 Note that the value of g NNr] is 0.07 in Table IV. 
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FIG. 7: (Color online) Same as Fig. 4 corresponding to the parameter set of Table IV. 



dictates to a large extent the behavior of the angular dis- 
tribution at forward angles. Figures 8(b,c) display the 
results corresponding to two additional parameter sets 
(not given here) with the same set of nucleon resonances 
as in Table IV. Overall, the two parameter sets yield 
a fit quality comparable to that of Fig. 7, but with very 



different resonance parameter values, which points to the 
ambiguity of the fit results if one relies solely on differ- 
ential cross-section data. Shown in Fig. 8(b) is the result 
obtained by using pure pseudovector coupling (A = 0) 
as in Fig. 7, while in Fig. 8(c) pure pseudoscalar cou- 
pling (A = 1) is adopted. The corresponding NNrj cou- 
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FIG. 8: (Color online) 77-meson angular distributions in the 
center-of-mass frame in "/p — > pr\ at T 7 = 2700 MeV and 
2887 MeV. (a) The results corresponding to the parameter 
set of Table IV, with QNN-q = 0.07. (b,c) The results obtained 
with two other parameter sets (not given in this work) using 
the pseudovector (A = 0) and pseudoscalar (A = 1) coupling 
choices, respectively, at the NNrj vertex. The resulting NNrj 
coupling constant values are QNN-q = 1-62 and 1.38, respec- 
tively. The contributions from the other nucleon resonances 
are practically negligible at these energies and are not dis- 
played. ./V15 is the sum of spin-5/2 resonance contributions. 
The data are from Crede et al. [1] . 



pling constants are g^^ v = 1-62 and 1.38, respectively. 
As can be seen in Figs. 8(b,c), in these cases the nucle- 
onic current contribution at backward angles is as large 
as that of the mesonic current at forward angles. How- 
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FIG. 9: (Color online) (a) Beam asymmetry Es and (b) 
target asymmetry Et in jp — > pr\ at T 7 = 2700 MeV and 
2887 MeV as a function of ^-emission angle in the center- 
of-mass frame. The solid curves represent the results ob- 
tained with the parameter set of Table IV, which correspond 
to Fig. 8(a). The dashed and dash-dotted curves are ob- 
tained with the two other parameter sets that correspond to 
Figs. 8(b) and (c), respectively. 



ever, the resonance contribution exhibits a similar 
angular dependence as the nucleonic current with a com- 
parable magnitude and interferes destructively with the 
nucleonic current. The destructive interference is almost 
complete and results in very small cross sections at back- 
ward angles as observed in the data. Overall, everything 
else being very similar between Figs. 8(b,c), we find no 
real sensitivity of the differential cross sections at these 
energies as to whether one employs pseudovector or pseu- 
doscalar couplings. 

The situation changes when one considers spin observ- 
ables. While there is no real difference at lower energies, 
however, at the high energies (T 7 = 2700 ~ 2887 MeV) 
considered here, the beam asymmetry Tib shown in 
Fig. 9(a) can distinguish clearly between the parame- 
ter sets corresponding to Fig. 8(a) on the one hand 
and Figs. 8(b,c) on the other. The marked differences 
are due to the marked differences in the values for the 
NNrj coupling constant g NNrj , which is vanishingly small 
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(g NNr , = 0.07) for the set corresponding to Fig. 8(a) 
and much larger (and about the same, g NNr) — 1-62 and 
g NNn — 1-38, respectively) for Figs. 8(b) and (c). This 
finding shows that the beam asymmetry at higher ener- 
gies can impose more stringent constraints, in particular, 
on the NNi] coupling constant. In any case, judging 
from the results in the present investigation, we expect 
the upper limit of the NNr/ coupling constant to be not 
much larger than g NNr) ~ 1.7. 

These parameter sets also lead to noticeable differences 
- albeit not quite as large — for the target asymme- 
try £t, as shown in Fig. 9(b). Of particular impor- 
tance, however, is that this observable can clearly dis- 
tinguish between whether one employs pseudoscalar or 
pseudovector coupling at the NNrj vertex which the re- 
sults of Figs. 8(b,c) and the respective curves of Fig. 9(a) 
cannot do. 

Before closing this subsection, we remark that, quite 
recently, the authors of Rcf. [11] have addressed the issue 
of chiral symmetry in ?7-meson photoproduction through 
the pseudoscalar-pseudovector mixing parameter A at the 
NNi] vertex by investigating this reaction close to the 
threshold. Our study reveals that one must be care- 
ful with such an investigation for the reasons mentioned 
above. In particular, if the NNrj coupling constant turns 
out to be very small, it will be very difficult to determine 
the value of the mixing parameter A. 



C. 7r p — > n,r\ 

In this subsection, we discuss the 7r~p — > nq reaction 
with the parameter sets determined above. The results 
for the total and differential cross sections for the re- 
action tt~p — ► nrj corresponding to the parameter sets 
of Tables I and II are displayed in Figs. 10(a,b) and 
Figs. 10(c,d), respectively. The cross-section results show 
that the various dynamical contributions of the two sets 
are very similar. The total cross section is rather well 
reproduced up to W rts 1.6 GcV, where it is dominated 
by the Su resonances, especially, by the Sii(1535) res- 
onance. Here, both the nucleonic and mesonic currents 
yield very small contributions. However, we do not re- 
produce the total cross section at higher energies due to 
the absence of the higher-mass resonances in these pa- 
rameter sets and the absence of the Nirir contribution 
via the coupled channel [30, 32] in this model. For differ- 
ential cross sections, we again note that these parameter 
sets are unable to reproduce the structure exhibited by 
the data at higher energies. 

As we have shown before, the inclusion of the 
Pi3(1720) does not improve the results for the jp — > pi] 
reaction significantly. However, this resonance provides 
an important contribution to reproduce the structure ex- 
hibited by the differential cross section data in n~p — > 
nq. This is illustrated in Fig. 11 corresponding to the 
parameter set of Table III. In addition, the Pi3(1720) 
resonance also helps improve to some extent the fit qual- 



ity for the total cross section above W « 1.6 GeV. 

Figure 12 shows an alternative fit for the total and dif- 
ferential cross section using the same set of resonances 
as those of Fig. 11 plus Pn(1710). The corresponding 
fit results for photoproduction are of comparable qual- 
ity to those shown in Sec. Ill A. As one can see from 
Fig. 12(a), the small bump near W = 1.7 GeV in the 
spin- 1/2 resonance contribution (dashed line) is caused 
by the Pn(1710) resonance which makes the bump in 
the total contribution (solid line) more pronounced com- 
pared to the result of Fig. 11(a). The Pu(1710) reso- 
nance seems also to affect the differential cross section 
in the vicinity of W = 1670 MeV, improving the agree- 
ment with the data to some extent. It is interesting to 
note that the recent chiral constituent quark-model cal- 
culation of Ref. [35] shows a dominant contribution from 
the Pn(1710) resonance around W = 1700 MeV, in con- 
trast to the present approach, where the dominant con- 
tribution arises from the spin-3/2 resonances. However, 
the authors of Ref. [35] have also found that the agree- 
ment with the measured differential cross sections in the 
W = 1609-1670 MeV energy range can be improved if 
the sign of their Pn(1710) partial wave amplitude is re- 
versed and if they employ a larger total decay width of 
T R ~ 350 MeV. The effect of the Pn(1710) resonance is, 
then, very similar to that found in the present calcula- 
tion. We note that the total decay width of the Pn (1710) 
resonance in Fig. 12 is r# = 95 MeV. 



D. NN -> NNrj 

In this subsection, we turn our attention to the 
NN — ► NNrj reaction. Although the model calculation 
of Ref. [26], which is based on a strong P-wave con- 
tribution, reproduces nicely the shape of the measured 
pp invariant mass distributions, it largely underestimates 
the total cross section data near the threshold. Here we 
present new results based on a combined analysis of the 
reactions listed in Eq. (1). 

Shown in Fig. 13 are the results for the NN — > NNi] 
reaction corresponding to the parameter set of Table I. 
It can be seen that the present model reproduces rather 
reasonably all the considered data for this reaction, in- 
cluding the energy dependence of the total cross sections 
at lower energies. This is a considerable improvement 
over the results of Ref. [26]. However, we still underes- 
timate the total cross section in pp collisions by a factor 
of <~ 1.5 when the excess energy Q is Q < 10 MeV. We 
expect that the inclusion of prj FSI and possibly the three- 
body effects would resolve this discrepancy once they are 
properly taken into account. Here, the major difference 
of our results from the previous calculation of Ref. [26] 
is that we have a much stronger spin-3/2 resonance con- 
tribution to the cross sections, especially at lower ener- 
gies. This is due to the large coupling of the P13 reso- 
nances to the p and u> vector mesons, resulting from the 
global fitting procedure (cf. Table I). In contrast, the 
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FIG. 10: (Color online) Results for tt p — > nrj corresponding to the parameter set of Table I. (a) Total cross section as a 
function of the total energy of the system W. (b) r\ angular distribution in the center-of-mass frame. Here, iVll stands for 
5n(1535) + Silliest)) contributions and iV13 = Di 3 (1520) + Di 3 (1700) contributions. (c,d) Same as (a) and (b) but with the 
parameter set of Table II. The numbers in (b) and (d) denote the total center-of-mass energy W in MeV. The labeling of the 
curves in (b,d) is the same to that of (a,c). The experimental data are from Ref. [27, 28]. 




FIG. 11: (Color online) Same as Fig. 10 but corresponding to the parameter set shown in Table III. Here, 7V13* = Z?i3(1520) + 
Dis ( 1700) +Pis( 1720). 
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Sii resonance contribution to the cross section is surpris- 
ingly small, especially at lower energies, which is due to 
the strong destructive interference among the exchanged 
mesons in the excitation of Sn resonances. Note that in 
the present calculation the resonance coupling constants 
involving vector mesons, g R ^y {V — p,uj), are basically 
fixed by the NN — > NNr) reaction, while the coupling 
constants to pseudoscalar mesons, g RNP [P — 7T, rf), are 
fixed by the jp — > prj and n~p — > nq reactions to a large 
extent. We also found that the mesonic current yields a 
very small contribution to the cross sections. The nuclc- 
onic current contribution is also negligible because of the 
very small NNr] coupling constant that results from the 
fit to the photoproduction data, as discussed above. 

It should be emphasized that it still remains to be ver- 
ified whether the dominance of the D13 resonances dis- 
cussed above is indeed true. In fact, in spite of the present 
lack of information on the corresponding coupling con- 
stants, g RNV , the obtained values (cf. Table I) may be 
too large to be realistic. For example, a rough estimate 
of these coupling constants from the PDG hclicity ampli- 
tudes [15] in conjunction with vector-meson dominance 
yields 



My, i.e., by writing 



for R = Di 3 (1520): 



and 



for R = D 13 { 1700): 
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where ^ RNV is the effective hadronic coupling function 
that includes the form factor F given by Eq. (A. 17). 
We emphasize in this context that these coupling con- 
stants cannot be determined uniquely in the present anal- 
ysis. Indeed, as shown in Fig. 14, a scenario in which 
the S\\ resonance dominates over the D13 resonance 
current can be achieved. Therefore, the consistency of 
these coupling constant values with other independent 
reaction processes should be examined in detail. For 
this purpose, vector-meson production processes, such 
as NN -> NNV, nN -> NV and jN -> NV, are of 
particular interest [50]. 

In Ref. [20], where the dominant 77-production mecha- 
nism is the excitation of the 5*11(1535) resonance via the 
pion-exchange, the analyzing power A y exhibits a zero 
at cos(^) ~ 0. In the present calculation, where the 
dominant production mechanism is the D13 resonance 
excitation, the zero of A y is shifted toward backward an- 
gles at larger Q. As has been pointed out in Ref. [20], 
unlike the total and differential cross sections, the ana- 
lyzing power is very sensitive to the reaction dynamics. 
In fact, the dashed curves in Fig. 13 for A y represent the 
results of another parameter set (not given here) which 
yields practically the same results for other observables 
considered here (see also the results shown in Fig. 14 be- 
low) . Unfortunately, the data are not accurate enough to 
disentangle these different scenarios. More accurate data 
will, therefore, impose more stringent constraints so as 
to help distinguish different dynamics of 77 production in 
NN collisions. 

The results corresponding to the parameter set of Ta- 
ble II are shown in Fig. 14. They are of comparable 



These values are corrected to the normalization point, at 
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FIG. 13: (Color online) Results for the reaction NN — > NNrj corresponding to the parameter set of Table I. (a) Total cross 
section as functions of the excess energy Q in pp and pn collisions, (b) i) angular distribution in the overall center-of-mass frame, 
(c) Final proton angular distribution, (d) pp invariant-mass distribution, (e) prj invariant-mass distribution, (f) Analyzing 
power. In (f), only the total contributions are shown (solid curves); the dashed curves represent the correspond results of 
another parameter set (not given in this work) which yields practically the same results for other observables considered in this 
reaction. The labeling of the curves in (b,c,d,e) is the same as that in (a). The data are from Refs. [16-19]. 



fit quality to those in Fig. 13 overall. However, the dy- 
namical content is quite different. In this case, the Sn 
resonance dominates over the -D13 resonance contribu- 
tions to the cross sections. Unlike the results shown in 
Fig. 13, there is no strong destructive interference among 
the exchanged meson contributions to excite the interme- 
diate S11 resonances. Of course, the smallness of the D13 
resonance contributions is directly correlated to the very 
small constants for the coupling of the D\% resonances 
to vector mesons (cf . Table II) . The analyzing power for 
this parameter set at Q — 40 MeV exhibits a qualita- 
tively different behavior from that found in Fig. 13. But 
again, the quality of the current experimental data does 
not allow to draw any definite conclusion. 

It should be emphasized that the much larger differ- 



ence in the dynamics in NN — > NNrj between the two 
parameter sets considered above as compared to those 
in 7p — > prj and 7r~p — > nr\ stems from the much richer 
interference effects in the former reaction. In particular, 
note that, for each meson M exchanged (jr, ?y, p, u), there 
are two coupling vertices involved in the nucleon reso- 
nance currents (i.e., RNM and NNM, see Fig. 3). This 
allows for an interference among all the exchange mesons 
involving the RNM coupling constants, a feature that is 
absent in the other two reactions. Therefore, the meson- 
production reactions in NN collisions, in conjunction 
with other basic photon- and (two-body) hadron-induced 
reactions, should help in extracting the resonance pa- 
rameters to a large extent. In particular, adding the 
vector meson (V) production channels, NN — > NNV, 
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FIG. 14: (Color online) Same as Fig. 13 for the parameter set of Table II. 



Q=74MeV / ■■ 




FIG. 15: (Color online) Predictions for the invariant-mass 
distributions da/drn^ x for pp (x — p) and pr\ (x = 77) in 
the reaction pp — > pprj corresponding to the excess energy of 
Q = 74 MeV for the parameter sets of Tables I (solid line) 
and II (dashed line). 



irN — » NV and jN — > NV, to the list of reactions given 
in Eq. (1), should impose more stringent constraints on 



the resonance coupling constants. 

Finally, Fig. 15 shows predictions using the parameter 
sets of Tables I and II for the invariant-mass distributions 
for pp and p-q at the excess energy of Q = 74 MeV. At 
this energy, data have been measured at CELSIUS with 
the WAS A detector [51]. 



IV. SUMMARY 

In the present work, a combined analysis of the re- 
actions 7p — > pr/, ir~p — > 77,77, PP ~^ PPVi an d W 1 ^ pnr) 
has been carried out within a relativistic meson-exchange 
model of hadronic interactions. Both jp — ► pr\ and 
n~p — > nq reactions have been treated in the tree-level 
approximation with the former reaction containing a gen- 
eralized contact current which ensures gauge invariancc 
of the reaction amplitude. The NN — > NN-q reaction 
has been treated in the DWBA approximation with the 
explicit treatment of the NN ISI and FSI. The free pa- 
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ramctcrs of the model, especially the nucleon resonance 
parameters, are then fixed in a global fitting procedure. 

Overall, the photoproduction data can be described 
reasonably well with the inclusion of the well-established 
511(1535), 511(1650), £> 13 (1520), and L>i 3 (1700) reso- 
nances as the minimally required set of resonances to 
achieve a reasonable fit to the currently available data. 
The inclusion of additional well-known resonances in the 
same mass region [including the spin-5/2 £>i5(1675) and 
Fi5(1680) resonances] does not further improve the qual- 
ity of the overall description of the data. The measured 
angular distributions at higher energies and backward 
angles are compatible with a vanishing NNrj coupling 
constant. However, in order to extract this coupling con- 
stant unambiguously within an approach of the type pur- 
sued here, one needs to go beyond the resonance region to 
avoid possible interference effects. On the other hand, as 
we have seen in the results of Fig. 9, the beam asymme- 
try can impose constraints on the NNrj coupling constant 
that are much more stringent than the differential cross 
sections because of the interference between the nucleonic 
and resonance currents. One difficulty of the present cal- 
culation is to reproduce the sin(20,j) dependence exhib- 
ited by the measured target asymmetry near threshold 
and this shows that further investigations are necessary 
to understand better the production mechanism of this 
reaction. 

Our model can also explain the available data on 
tt~P — > nrj reasonably well, at least for energies not too 
far from the threshold. At higher energies the total cross 
section is underestimated due to the lack of the Nttt: 
contribution via the coupled channel. 

Our model also describes the NN — ► NNrj reaction 
data rather well. This is in contrast to the results of 
Rcf. [26], where the pp and pij invariant mass distribu- 
tions are well-described but the ppr\ total cross section 
is underestimated by a factor as large as about 5 near 
threshold. This problem has been cured to a large extent 
in the present approach. However, we still underestimate 
the total cross section by a factor of ~ 1.5 in the region of 
Q < 10 McV, which is expected to be addressed once the 
Ni] FSI and three-body effects are properly taken into 
account. We also verified that the analyzing power is 
sensitive to the reaction dynamics. Unfortunately, how- 
ever, the currently available data are not accurate enough 
to unambiguously distinguish different dynamical contri- 
butions. 

As we have illustrated with some selected examples, 
the present approach is unable to determine a unique 
set of (resonance) parameter values. In fact, it was 
shown that different parameter sets which describe the 
data equally well lead to results which exhibit quite dif- 
ferent reaction dynamics. This is especially true for 
the NN — > NNr/ reaction, where we found quite dif- 
ferent values of the RNV coupling constants. The in- 
clusion of vector-meson production reactions into the 
combined analysis should help constrain those coupling 
constants. Consequently, our study reveals that one 



must be cautious in interpreting the resonance param- 
eters extracted from this kind of analyses, especially, if 
only a single reaction process is considered. It is clear 
that in order to extract more accurate information on 
the nucleon resonances, one must combine the investiga- 
tion of hadron-induced meson production with the cor- 
responding photon-induced reactions. To help in this, 
ample data sets are now available for photo- and electro- 
production processes from various accelerator facilities 
including the Thomas Jefferson National Accelerator Fa- 
cility, SPring-8, CB/ELSA, GRAAL, etc. By contrast, 
data for hadron-induced production processes are much 
more limited and we clearly need more of them. 

In this connection, we note that the progress in the 
study of meson-production processes in NN collisions, 
both experimentally and theoretically, has reached such a 
level that it allows us to address certain concrete physics 
issues, especially, when they are investigated in conjunc- 
tion with other independent reactions. This has been 
illustrated in the present work for the specific case of 
r] production where some information on nucleon reso- 
nances can be extracted. In particular, the consideration 
of meson-production processes in NN collisions in con- 
junction with photon- and (two-body) hadron-induced 
reactions aiming at a resonance parameters extraction, 
should help impose more stringent constraints on these 
parameters. As pointed out in the last subsection, me- 
son production in NN collisions exhibit much richer in- 
terference effects, a feature that is absent in more basic 
two-body reactions. Furthermore, the inclusion of this 
reaction in the resonance-parameter extraction is espe- 
cially relevant because the existing data for meson pro- 
duction (other than for the pion) in two-body hadronic 
reactions are rather scarce and of relatively low accuracy. 
Currently, there exist only very limited efforts to improve 
or extend the corresponding database [28] . On the other 
hand, the available data on meson production in NN 
collisions arc much more accurate and, moreover, the 
corresponding database can be and is being expanded, 
especially at the COSY facility. 
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APPENDIX 

In this Appendix we give the ingredients which de- 
fine our models described in Sec. II. Throughout this 
paper, we use the notation N and R for the nucleon 



20 



and nucleon resonance fields, respectively; Mb denotes 
the mass of the baryon B (= N,R). We also use S 
(= a,a ), P (= 77, 7f), and (= uj^,p^) to denote 
the scalar, pseudoscalar, and vector meson fields, re- 
spectively. The vector notation refers to the isospin 
space. For isovector mesons, S = S ■ f, P = P ■ f, and 
Vn = • t. The mass of the meson M is denoted by Mm 
(M = S,P,V). The photon field is denoted by A^. We 
define V» v = d»V v - d v V» and = d»A v - d v A». 

We use the superscript j p in the Lagrangian densi- 
ties ') involving the nucleon resonance R to denote 
the spin-parity j p of that resonance. Furthermore, for 
convenience, wc define 



nucleon interactions are 



£nns = 9nns NNS, 
£nnp = -9nnpN |r (+) 



iX- 



£nnv — —Qnnv^ 



{ 



1 - A 
" 2M N 



2M 



N 



(A.2) 
P}N, (A.3) 



V, N, (A.4) 



r (+) 
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and 



r(-) 



(A.l) 



1. Hadronic interaction Lagrangians 

where the parameter A was introduced in Ln np to inter- 
The following interaction Lagrangian densities describe polate between the pseudovector (A = 0) and the pscu- 
thc hadronic vertices. The Lagrangians for meson- doscalar (A = 1) couplings. 

I 



The effective Lagrangians describing the interactions of nucleon resonance with the nucleon and pseudoscalar meson 
P or vector meson V read 



£rnp — T9rnpR \ r (±) 



'-RNV — ~ - - 1XL 



2M N 



iX + 



9rnv 



1 - A 



M R ±M N 

( 7^ 2 
\MrtM n 



PjN+ H.c., 



V>*}N + H.c, 



(A.5) 
(A.6) 



for a resonance of spin-|. The Lagrangian (A.5) contains the pseudoscalar-pseudovector mixing parameter A, similar 
to Eq. (A.3). 

For spin- 1 resonances, we use 



^RNP = g - M yR IX ®A^ (± \d v P)N+ H.c. 



(A.7) 



,(2) 



(3) 



c ( lP v = -i^R f3 e^ lu v^N - Vj^fRPe^v^dvN T fp^e^r^ {d u vn n + h. c, (a. 



where the coupling tensor is 6 M1 , = g^ v — (z + 5)7^7,, and we take the off-shell parameter z = — \ for simplicity. 
For the interaction Lagrangians of spin-| resonances, we use 



^np = ^-^f^ {±) i^d v P) R, v + H. c, 



(A.9) 



,(2) 



,(3) 



+ H.c, 



pa 

(A.10) 



The hadronic interaction Lagrangians among mesons 
are 



= -^s a ^(d a p") ■ {d v p»)r, , (A.lla) 



e a ^{d a ^){d»uj») n , (A.llb) 



V y/M v M v ^ 



(A.llc) 
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with the convention £0123 = — 1 for the Levi-Civita anti- 
symmetric tensor. 

We follow Refs. [20, 38] for the coupling constant values 
at the NNM vertices above, except for the coupling con- 
stant g NNn , which is treated as a free parameter in the 
present work. The coupling constants of the RNM inter- 
action Lagrangians, as well as the resonance masses Mr, 
are free parameters to be adjusted to reproduce the data. 
For the other coupling constants, following Ref. [20], we 



use g. 



■qpp 



4.94, 



4.84, and g r 



1.81. 



interaction of the nucleon reads 



-AW 7 



= -eN 



ey 



n _ 



a» v d„ 



2M 



N 



N. 



(A.12) 



2. Electromagnetic interaction Lagrangians 

The electromagnetic vertices are calculated from the 
Lagrangian densities given below. The electromagnetic 



where e stands for the elementary charge unit, and e = 
(1 + t 3 )/2 and k = K p (l + t 3 )/2 + K n (l- r 3 )/2, with the 
anomalous magnetic moments n p — 1.739 of the proton 
and k„ — —1.931 of the neutron. 



The photo-transition Lagrangians of resonances into the nucleon are 

= e^RT^a^A^N + H. c, (A.13) 

4n] = -ie^R^T^F^N - e^Bfe^F^N + H. c, (A.14) 

(1) (2) 

= tmpyi ^> rW ( dap n R ^ ^mPv^)^ ( dap ^ R ^ + H - c " (A.15) 



and the photo-transition Lagrangian between mesons are 
Cnry - -e^e^^pgWAnv, (A.16a) 



C 



„,n = -e^e^id^WA^. (A. 16b) 



The electromagnetic coupling constants involving nu- 
cleon resonances in the above Lagrangians are free pa- 
rameters to be adjusted to fit the pertinent experimen- 
tal data. The coupling constants g^ vl (v = p, lu) in the 
above equations are determined from a systematic anal- 
ysis of the pseudoscalar and vector meson radiative de- 
cays [20, 52], which leads to g vpp = 1.44 and 3 = 0.47. 
Their signs are inferred also from SU (3) flavor symme- 
try considerations in conjunction with the sign of the 
coupling constant g vvl determined from a study of pion 
photoproduction in the 1 GeV energy region [53]. 



Form Factors 



where p' and p denote the four-momenta of the two 
baryons, and q is the four-momentum of the meson at 
the three-point vertex. Here, we use 



F B (x) 



K% + {x-M%f 



(A.18) 



where the cutoff A# = 1200 McV is taken the same for 
all the baryons, and Fm{q 2 ) is given by 



M{q ' ~ A 2 -„2 



A 2 _ M 2 \ « 
M 1V1 M 



(A.19) 



with n — 1 for a scalar or a pseudoscalar meson and 
n = 2 for a vector meson. The values of Am are taken 
the same as those used in Ref. [20]. 

The electromagnetic rjvj vertex is multiplied by the 
form factor G v (q 2 ) which describes the off-shell behavior 
of the intermediate vector meson with squared momen- 
tum transfer q 2 (cf. the fourth diagram in Fig. 1). In 
general, we use the dipole form 



Each vertex obtained from the interaction Lagrangians 
given above is multiplied by a phcnomenological cutoff 
function 

F(p' 2 ,p 2 ,q 2 ) = F B (p' 2 )F B (p 2 )F M (q 2 ), (A.17) 



G v (q 2 ) = 



A 



*2 



Af-q 2 ) ' 



(A.20) 



where the cutoff A*, taken to be same for both p and w, 
is a free parameter. 
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4. Propagators 

The calculation of the Feynman diagrams displayed in 
Figs. 1, 2, and 3 requires the corresponding baryon and 
meson propagators. Their explicit forms are given here. 



S(p) 
A(<7 2 ) 
D^(q) 



1 



M N 
1 



q 2 - M 2 M + le ' 
-g» v + q^q u /M 2 



q 2 - Mi 



(A.21) 
(A.22) 

(A.23) 



where S(p) is the nucleon propagator with the nuclcon 
four-momentum p, A(q 2 ) is the scalar or pseudoscalar 
meson propagator with four-momentum q, and D^ u (q) is 
the vector meson propagator. 

For a spin-1/2 resonance propagator, we use the ansatz 



Sy 2 (p) 



1 



$ - Mr + § 



r 

Mr 



p 2 -M 2 R + Ui> + M R )T 1 



(A.24) 



where T is the energy dependent resonance width. 

For spin-3/2, the resonant propagator reads in a 
schematic matrix notation 



Similarly, the propagator for a spin-5/2 resonance is 
given by 



#5/2 0) 



G. 



(A.27) 



where 



= G /3a 



i 5 



8f3o 



with 



(A.28) 



p^j> 
Ml' 



(A.29) 



The above spin-5/2 propagator is a variant of the one 
given in Ref. [54]. 

We write the resonance width T as a function of W — 
i/i according to 



T{W) = T R 



N W 7 

i=i j=i 



, (A.30) 



Sz/2<J>) 



(j/> - M R )g 



A ' 

? 2 r 



A , 



(A.25) 



where all indices are suppressed, i.e., g is the metric ten- 
sor and A is the Rarita-Schwinger tensor written in full 
detail as 

A ^ = " 5 ^+3^e7e. + ^^ + ^ , 

(A.26) 

where a, (3, and e enumerate the four indices of the 7- 
matrix components (summation over e is implied). The 
inversion in Eq. (A.25) is to be understood on the full 
16-dimcnsional space of the four Lorentz indices and the 
four components of the gamma matrices. The motivation 
for the ansatz (A.25) and the technical details how to 
perform this inversion is given in Ref. [38]. 



where the sums over i and j respectively account for de- 
cays of the resonance into N and two- or three-hadron 
channels and into iV 7 radiative decay channels. The total 
static resonance width is denoted by Tr and the numer- 
ical factors f3 i and 7- (with < (3 i ,^j < 1) describe the 
branching ratios into the corresponding decay channel, 
i.e., 



N N-, 



(A.31) 



i=l 



We parametrize the width functions Ti and T 7 . (which 
are both normalized to 1 at W = Mr) as given in 
Ref. [38] to provide the correct respective threshold be- 
havior. 
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